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Abstract 

Hyperfine intervals in light hydrogenic atoms and ions are among the most 
accurately measured quantities in physics. The theory of QED corrections has 
recently advanced to the point that uncalculated terms for hydrogenic atoms 
and ions are probably smaller than 0.1 parts per million (ppm), and the ex- 
periments are even more accurate. The difference of the experiments and QED 
theory is interpreted as the effect on the hyperfine interaction of the (finite) 
nuclear charge and magnetization distributions, and this difference varies from 
tens to hundreds of ppm. We have calculated the dominant component of the 
Is hyperfine interval for deuterium, tritium and singly ionized helium, using 
modern second-generation potentials to compute the nuclear component of the 
hyperfine splitting for the deuteron and the trinucleon systems. The calculated 
nuclear corrections are within 3% of the experimental values for deuterium and 
tritium, but are about 20% discrepant for singly ionized helium. The nuclear 
corrections for the trinucleon systems can be qualitatively understood by in- 
voking SU(4) symmetry. 
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1 Introduction 



The physics of hyperfine structure (hfs) is driven by magnetic interactions. This 
physics has a short-range nature, and is more complicated and challenging than 
"softer" regimes in atomic physics. This is especially true of the nuclear contribution 
to hyperfine structure, because the nuclear current density is more complicated than 
the nuclear charge density and less well understood [1, 2]. 

Until very recently hyperfine splittings in light hydrogenic atoms were by far the 
most precisely measured atomic transitions. Many of these very accurate experiments 
date back nearly half a century. Theoretical predictions are far less accurate, but have 
improved considerably in recent years. Non- recoil and non-nuclear contributions [3, 4] 
are known through order a 3 Ep, where Ep is the Fermi hyperfine energy (viz., the 
leading-order contribution) and a is the fine-structure constant. Because the hadronic 
scales for recoil and certain types of nuclear corrections are the same, recoil corrections 
are treated on the same footing as nuclear corrections [3], and we will call both types 
"nuclear corrections." It is very likely that the uncalculated QED terms of order 
a 4 Ep in light atoms contribute less than 0.1 ppm. Although much larger than the 
experimental errors, this is still significantly smaller than the nuclear corrections. We 
restrict ourselves to hydrogenic s-states in this work, because these states maximize 
nuclear effects. 

Table I. 

Table 1: Difference between hyperfine experiments and QED hyperfine calculations 
for the nth s-state of light hydrogenic atoms times n 3 , expressed as parts per million 
of the Fermi energy. This difference is interpreted as nuclear contributions to the hy- 
perfine splitting [4] . A negative entry indicates that the theoretical prediction without 
nuclear corrections is too large 



n 3 ( £exp _ E Q™ )/Ef (ppm) 



State 


H 


2 H 


3 H 


3 He+ 


Is 


-33 


138 


-38 


-212 


2s 


-33 


137 




-211 



Table I is an updated version of the corresponding table in Ref. [4] . Because hyper- 
fine splittings are dominated by short-range physics, we expect the splittings in the 
nth s-state to be proportional to |0 n (O)| 2 ~ 1/n 3 , where <p n {r) is the non-relativistic 
wave function of the electron. Forming the fractional differences (in parts per million) 
between and E^ D leads to the tabulated results. These large differences reflect 
neither experimental errors nor uncertainties in the QED calculations; they directly 
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reflect large nuclear contributions. 

One complication in performing the nuclear calculation is obtaining a final re- 
sult that is tractable for numerical calculations. A framework fortunately exists for 
performing systematic expansions[2] of nuclear matrix elements in powers of (Q/A), 
where Q is a typical nuclear momentum scale that can be taken to be roughly the 
pion mass (m n ~ 140 MeV), and A is the large-mass QCD scale (~ 1 GeV) typical 
of QCD bound states such as the nucleon, heavy mesons, nucleon resonances, etc. 
This framework, called power counting, also extends to nuclei, where 1/Q specifies a 
typical correlation length (and a reasonable nearest-neighbor distance) in light nuclei 
(~ 1.4 fm). This expansion in powers of (Q/A ~ 0.1 — 0.15) should converge moder- 
ately well. For the purposes of this initial work, we restrict ourselves to leading-order 
terms in the nuclear corrections. 

This restriction eliminates nuclear corrections of relativistic order, which we will 
briefly discuss later. There have been relatively few relativistic calculations in light 
nuclei because of the complexity of the nuclear force. The calculations that exist are 
known to generate rather small corrections [5], and they are especially small in the 
deuteron case because of its weak binding[6]. Few calculations exist for the much 
more complicated three-nucleon systems [7]. Most nuclear physics knowledge and lore 
in light nuclei is non-relativistic in nature. 

In processes that involve virtual excitation of intermediate nuclear states (each 
state \N) with its own energy, E N , relative to the ground-state energy, E ) the ex- 
citation energy (un = E N — E Q ) is of order Q 2 / A and typically is a correction to 
the leading order [2]. Consistency therefore demands that we drop such terms. The 
nuclear recoil energy Q 2 /2M, where M is the nuclear mass, has the same scale and 
can also be dropped. 

Our goal is to evaluate the deuteron contribution in leading order as carefully as 
possible, and to use the impulse approximation to evaluate the 3 He and 3 H results. 
This restriction should be accurate to within the 10-15% uncertainty of our leading- 
order approximation and adequate for a first calculation of the trinucleon sector. This 
level of accuracy does require the inclusion of the intrinsic hyperfine structure of the 
nucleons. We will see that the calculated results are good for deuterium and tritium 
and fairly good for 3 He + . Our goal is to present a simple and compelling picture of the 
nuclear hyperfine structure in d, 3 H, and 3 He + , using a unified approach for all. To 
accomplish this we invoke a simple and intuitive model of trinucleon structure based 
on SU(4) symmetry that is sufficiently accurate to explain the patterns in Table I. 
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2 Nuclear Contributions to Hyperfine Structure 



The hyperfine interactions that interest us are simple (effective) couplings of the elec- 
tron spin to the nuclear (ground-state) spin: a • S where a is the electron (Pauli) spin 
operator and S is the nuclear spin (i.e., the total angular momentum [8]) operator. 
Other couplings of the electron spin are possible and either generate no hyperfine 
splitting, none in s-states, or higher-order (in a) contributions. Coupling of the elec- 
tron spin to the electron angular momentum (i.e., the electron spin-orbit interaction, 
which vanishes in s-states) is one example. 

We begin with a sketch of the Fermi hyperfine splitting in order to establish our 
notation and conventions. The electron charge operator (ip'ijj) is a space scalar, and 
by itself does not generate a term of the desired form because nuclear-spin informa- 
tion cannot be transmitted via the exchange of the appropriate component of the 
virtual photon's propagator. The leading-order term in the hyperfine energy shift 
consequently has the (usual) form of a current-current interaction. 

E F = a J d 6 r J d 3 r' r v y r v ^ — , (1) 

where ip^aip is the current density of the (Dirac) electron and J(r') is the nuclear 
(ground-state) current density. Ignoring higher-order (in a) terms this expression can 
be manipulated into the form 



2m, 



[ d 3 r |0 n (r)| 2 / dh-'J(r') • a x V r — , (2) 
J J r — r 



where m e is the electron mass and (f) n (r) is the usual non-relativistic atomic nth 
s-state wave function. The restriction to atomic s-states limits the current to the 
magnetic-dipole part, which can be written in terms of the magnetic-dipole operator, 
jl, as 

A r) -> J M i( r) = V r x /2( f) , (3a) 

and jl can be written in terms of the nuclear ground-state magnetization density 

p M (r) as 

— * 

it(f) = ^PM(r), (36) 

where / d 3 r' Pm{t') = 1 and (SS\fi z \SS) = fi N defines the nuclear magnetic moment. 
We finally obtain the Fermi hyperfine energy in the limit Pm{^') 5 3 (r'): 

E F = ^ a -A j Sr\Ur)?PM(r) = ^'^f^ ^ ■ (4) 
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where the factor of (a ■ S)/S leads to a hyperfine splitting of (2 S + 1)/S. All ad- 
ditional contributions will be measured as a fraction of this energy. We note that 
using the point-like pu and Eqn. (3a) in Eqn. (1) leads to an atomic matrix element 
/ °° drF(r)G(r) (in the standard Dirac notation) that can be evaluated analytically 
for higher-order Coulomb corrections [3] to the leading contribution of 0(Z 3 a 4 ), where 
Z is the nuclear charge. Note that E F is independent of nuclear structure. 





(a) 

Figure 1: Nuclear Compton amplitude with direct (a), crossed (b), and seagull (c) 
contributions illustrated. Single lines represent an electron, double lines a nucleus, 
wiggly lines a photon propagator (with four-momentum </*), and shaded double lines 
depict a nuclear Green's function containing a sum over nuclear states. The seagull 
vertex maintains gauge invariance. 

Naively calculating higher-order (in a) corrections from the first part of Eqn. (4) 
fails because only that part of (f> n (r) inside the nucleus (i.e., within the magnetic 
"size", R M ) contributes to the integral, and that is the (only) part of the electron's 
wave function significantly modified by the nuclear charge distribution, p c h(r), whose 
radius, R c h, is nearly the same as Rm- In other words a proper calculation [9, 10] 
must take into account modifications of 0„(r) by p c h, and that necessarily involves 
one order in a higher than E F . Thus we need to perform a consistent second-order 
(in a) calculation of the electron-nucleus interaction, depicted in Fig. (1). 

These graphs are constructed from the direct, crossed, and nuclear seagull contri- 
butions to the nuclear Compton amplitude. Only the forward-scattering part of this 
amplitude is required for the 0(a) corrections to E F , and this generates a short-range 
atomic operator that samples the (upper-component) s-state wave functions only near 
the origin. The resulting energy shift is then given by 

A£ = i(4™)U(0)| 2 f ^- , - ^# f ^ ~ A w (5) 
V ; ir " V n J {2ny{q 2 + ie) 2 {q 2 -2m e q + ie)' V ; 

where is the lepton Compton amplitude and T^ v is the corresponding nuclear 
Compton tensor, both of which are required to be gauge invariant. The lepton tensor 
can be decomposed into an irreducible spinor basis; we can ignore odd matrices and 
spin-independent components, since they do not contribute to the hyperfine structure 
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at 0(aEp). We also ignore (with one exception treated later) terms that couple two 
currents together. It is easy to show that since the nuclear current scales as 1/A (the 
conventional components of the current have explicit factors of 1/M), two of them 
should scale as 1/A 2 and generate higher-order (in 1/A) terms. This leaves a single 
term representing a charge-current correlation 

AE - (4na) |0 n (O)| J — ^ - - — - - } . (6) 

The seagull terms B 0m (q,—q) and B m0 (q,—q) are of relativistic order [11] (~ 1/A 2 ) 
and can be dropped. Although the term B mn (q, — q) is of non-relativistic order, it 
does not contribute in conventional approaches because of crossing symmetry. 

The explicit form for T m0 (suppressing the nuclear ground-state expectation value, 
but including all intermediate states, N) is 

T m0 (g, -q) = V (ZtiMM + P(v)\ N )( N \ jm (-v) ) ) (?) 
n V q -^N + ie -go - u N + ie J 

which greatly simplifies Eqn. (6) in the limit ojn — > and m e — > 0, leading finally in 
this limit to a very simple form when closure is used 

which is infrared divergent. Using 

J m (-q) = Jd 3 yJ m (y)e-^y, (9a) 

and 



p(q)= J d 6 x p(x)e iq - x , (9b) 
together with z = x — y, and a lower-limit (infrared) q— cutoff, e, we find 

AE = -8a 2 |0 n (O)| 2 / d 3 x J d 3 y{p(x),a-J(y)} x V z ^ - ^ + ^ , (10) 

where there is an implied (nuclear) expectation value. The constant term does not 
contribute because of the derivative, the z 2 -term contributes a term proportional to 
> anc ^ the ^ as ^ term is the term we are seeking. 

The remaining singular (second) term must be treated more carefully. The part 
proportional to y leads (because of the integral over x) to a contribution proportional 
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to Z, the total nuclear charge, and fi N , the nuclear magnetic moment. Since this 
contribution is already part of E F , keeping this term would amount to double count- 
ing, and we therefore ignore it. The x-term on the other hand generates (unretarded) 
dipole transitions, and the singularity (1/e) arises from neglecting 1/ujn and m e terms. 
Siegert's theorem[l] for unretarded electric dipole transitions (/ d 3 yJ(y) ~ iuj^D) 
generates an additional factor of u>n (i.e., e) so this singular term is actually of the 
form (e/e) — > (0/0) and requires a careful calculation. A rather tedious evaluation 
of this term in Eqn. (6) leads to AE a = -2ia 2 \(j) n (0)\ 2 a ■ D (In 2{H ~ m + 3/2) x D 
where H — Eq = ujn in intermediate states. A similar term arises in the product 
of two currents mentioned above Eqn. (6) and leads to AE^ = — ia 2 \<fi n (0)\ 2 <j ■ D 
(l n 2l^_l/2) Xj D. 

These contributions are completely unimportant, as we now demonstrate. The 

— * — * 

constant terms (including ln(m e )) are proportional to a ■ D x D, which vanishes for 
non-relativistic dipole operators because they commute. Replacing ujn by a constant 
(viz., the closure approximation) similarly vanishes. It is straightforward to show 
that the nuclear matrix element also vanishes in zero-range approximation, where one 
neglects the deuteron d-state and the potential in intermediate states (the structure 
of the dipole operator weights the tails of the wave functions, and this minimizes 
the effect of the intermediate-state nuclear potential). In perturbation theory it is 
possible to show that only the spin-orbit combination of potentials in intermediate 
states contributes (i.e., central and tensor terms cancel), and this small potential is 
of relativistic order (1/A 2 ), which we have previously agreed to ignore. One can also 
show that the non- vanishing deuteron contribution is proportional to f] 2 , the square of 
the d- to s-wave asymptotic normalization constant (~ (0.025) 2 ), which is extremely 
small. A recent brute-force numerical calculation confirms these estimates [12]. The 
terms AE a and AEb are therefore numerically negligible and can be ignored. 

Our final result in leading order is a relatively simple expression originally devel- 
oped in a limiting case by Low[13] for the deuteron, as sketched by Bohr[14] for the 
same system: 

AE^ = 7m2|0 2 " ( ° )|2 / d 3 x J d 3 y {p(x),a • z x J(y)} z + --., (11a) 

where both an atomic and nuclear expectation value is implied, but has been ignored 
in Eqn. (11a) and subsequent equations for reasons of simplicity. 

A more convenient representation of this result is obtained by dividing both sides 
by the expression for the Fermi hyperfine energy given by Eqn. (4). Since the Wigner- 
Eckart Theorem guarantees that the resulting form of Eqn. (11a) must be proportional 
to a ■ S/S (which cancels in the ratio), we arrive at a simple but powerful expression 
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for the leading-order contribution: 



AE^ = -2m e a5 Low , (lib) 

where 

5low = — [d 3 x [d 3 y{p(x),(zxj(y)) z }\x-y\, (11c) 

16/% J J 

and an expectation value is required of the z (or "3") component of the vector zx J(y) 
in the nuclear state with maximum azimuthal spin (i.e. S z = S). The intrinsic size 
of the nuclear corrections is given by (—2m e aR) = —38 ppm [i?/fm], where [i?/fm] 
is the value of the Low moment in Eqn. (11c) in units of fm. The results of Table I 
therefore suggest (correctly) that Low moments are on the order of a few fermi in 
light nuclei, which is quite sensible. 



3 Nuclear Matrix Elements 



We predicate our discussion with the deuteron in mind. Other nuclei can and will 
be treated mutatis mutandis. The isospin of the deuteron (T = 0) makes it a useful 
first case. We note that the nuclear physics in Eqn. (11) involves the correlation 
between the nuclear charge operator, p(x), and the nuclear current operator, J(y). 
If one inserts a complete set of states between these operators, there will be both 
elastic contributions (i.e., ground-state expectation values) that are called Zemach 
corrections[10], and inelastic contributions (called nuclear polarization corrections). 
Although we will calculate (or estimate) both types, it is much easier to calculate the 
sum of the two. 

The nuclear charge operator contains both isoscalar and isovector pieces, and is 
non-relativistic in leading order. We ignore relativistic corrections, as we discussed 
earlier. We therefore write the charge operator in the form 

p{x) = Y.e l {\x-x> l \), (12a) 
i=i 

where 

e i (\x-x^\)^p i pl h (\x-^)+h i p n ch (\x-x' j }) : (126) 

and 

Pi = ( > ( 12c ) 
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and 

(12d) 

This decomposes the ith nucleoli's charge operator into proton plus neutron parts. 
The densities p p ch and p\ are the intrinsic charge densities of the proton and neutron, 
respectively, while pi and hi are the proton and neutron isotopic projection operators, 
respectively. The coordinate x\ is the distance of the ith nucleon from the nuclear 
center of mass. We expect that the neutron charge density should play a very minor 
role, and we will find (later) that its contribution is only a few percent of that of the 
proton. Rather sophisticated models exist for the Fourier transform (i.e., the form 
factor) of f? ch [l7\. ' 

The nuclear current operator is more complicated, even if we ignore relativistic 
corrections (which we will). The problem is the mechanism underlying the nuclear 
force (viz., the exchange of charged mesons), which can also contribute to the nuclear 
current in the form of meson-exchange currents (MEC). These currents largely vanish 
for isoscalar transitions (such as the deuteron ground state) because there is no net 
flow of charge, but they can be sizable (10% - 20%) for isovector transitions. One can 
show that their contribution to the deuteron in Eqn. (11) almost entirely vanishes, 
and we will henceforth ignore these currents below Eqn. (13). We formally expand 
the current into convection, spin-magnetization, and meson-exchange parts 

J(V) = Jc(y) + V„ x p(y) + J MEC (y) + • • • , (13a) 

where 

&(fl = E{^>W-3l)} ( 13fo ) 

is the nuclear convection current, ffj is the ith nucleon's momentum in the nuclear 
center-of-mass frame, 

A 

tt(y) = J2°it i i(\y-Xi\) ( 13c ) 
i=i 

— * 

is the impulse-approximation magnetic-moment operator, Jmec is the nuclear meson- 
exchange current, and 

M\y-x'i\) =Pi^ P P P M(\y-x' i \)+hifi n p n M (\y-x' i \) (13d) 

is the nucleon magnetization density for the ith nucleon expressed in terms of pro- 
tons and neutrons separately. The quantities p p and p n are the (total) proton and 
neutron magnetic moments, while p p M and p%i are the intrinsic proton and neutron 
magnetization densities (normalized to 1). 
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Using Eqns. (12) and (13) the energy shift in Eqns. (11) can be evaluated. Rather 
than split Eqn. (11c) into Zemach terms (by inserting intermediate ground states 

— # 

between p and J) and polarization terms (by inserting intermediate excited nuclear 

— * 

states between p and J) , we will use the fact that our charge and current operators are 
each given by a sum over single-nucleon operators. Thus their product can be decom- 
posed into single-nucleon plus two-nucleon operators. These forms are particularly 
convenient to evaluate. We first write in an obvious notation that 



^Low - <$ S pj n + <5 L ow • (I 4 ) 



Note that the quantity 5^ in was not part of Low's original work, nor was there any 



'spm 



evidence at that time that such a term might be significant. We next use Eqns. (13) to 
manipulate the magnetization part of the current in Eqn. (11c) into the form (recall 
that z = x — y) 



C = — E I d * x [ d3 y \*-y\ x'SfrjiW - 41) - l( 3ff 



3 - zz- a jr 



(15) 

while the convection current can be reduced to 

Cw = -— ^— E / d 3 x [ d*y\x-mi**i)*M\z-Zi\)tA\y-m- ( 16 ) 

16p N M i d J J 

The one-body (i = j) part of the convection-current contribution vanishes upon inte- 
gration, as does the second (tensor) term in the magnetization contribution. Shifting 
the variables x and y by x\ leads to 

C = E ((r^Pi + (r)^n t ) a\ , (17) 
i=i \ p N p N I 



where 
and 



(r)^ =JdPxJ d*yp p ch (x)p p M (y)\x-y\ = 1.086(12) fm, (18a) 

(r)^ = Jd 3 xf d s yp n ch (x)p n M (y)\x-y\ (186) 

determine the proton and neutron parts of the one-body current. Note that the 
quantities (r)^ and (r)™ 2 ™ are the usual proton and neutron Zemach terms, and we 
have listed in Eqn. (18a) the value of the proton Zemach moment recently deter- 
mined directly from the electron-scattering data for the proton[15] (the neutron has 
not been evaluated). In numerical work described below we will use simple forms 
for the neutron and proton form factors: a dipole form for the proton charge and 
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magnetic form factors and the neutron magnetic form factor (Fo(g 2 ) = jrj^tjpzyz ) 
and a modified Galster[16] form for the neutron charge form factor (F G = jp^-j^s). 
To incorporate into our calculations the numerical value given by Eqn. (18a) we use 
j3 = 4.0285 fm -1 , which reproduces this value for the dipole case (see Appendix A for 
moments and correlation functions determined by this choice of form factors). The 
r.m.s. radius determined by this j3 is 0.86 fm, which is slightly smaller than the proton 
charge radius[17] but slightly larger than its magnetic radius[18] and thus represents 
an average value. The much smaller neutron moment, (r)j^ (see Appendix A), can 
be adequately represented using this value of (3 and A = 0.0190 fm 2 , which determines 
the neutron charge radius[19]. These numbers lead to (r)?£ = —0.042 fm, which we 
will use below. Because this value is such a small fraction of the proton result, the 
uncertainty in the neutron value plays no significant role. 

Equation (17) is still a nuclear operator, and its expectation value depends on 
the nucleus. We begin with the deuteron (which has T = 0) and this eliminates 
the T3 terms in pi and h{. The spin terms then sum to (o\ + (T 2 ) z , which is not the 
total angular momentum (it lacks L, the orbital angular momentum contribution). 
The expectation value of (<?i + <7 2 ) z is 2S' Z (1 — §-Pd), where S z is the z-component 
of the deuteron total angular momentum operator and P D is the amount of D-wave 
in the deuteron wave function (typically, slightly in excess of 5.6%). In the state of 
maximum S z this leads to 

= (l-lPn)( ! ±(r)% ) + !±(r)®) , (19) 

which is completely dominated by the proton. Note that the D-w&ve prefers to anti- 
align with the spin, which leads to the reduction in Eqn. (19). 

A similar (though more complicated) analysis is possible for 3 He and 3 H (see Ap- 
pendix B). The traditional (and very useful) decomposition of the trinucleon wave 
function uses representations of spin-isospin symmetry (viz., SU(4)). In addition to 
the somewhat larger D-state component (Pd ~ 9%), the significant .S-wave compo- 
nent comes in two distinct types: the dominant S-state (P$ ~ 90%) with a completely 
antisymmetric spin-isospin wave function and completely symmetric space wave func- 
tion, and the mixed-symmetry S'-state (P s > ~ 1%). The representations of SU(4) were 
used long ago to decompose contributions to the trinucleon magnetic moments, and 
this leads to[20, 21] 

- \pd)) , (20a) 



E 



l±Tf 



= 2S 



7 



l + r 3 



+ 



7 
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where 

7 = ^(P s ,-P D )^-0.06, (21a) 

£ = 1-^-^ = 0.86, (216) 

specify the two sign cases in Eqns. (20a) and (20b). Note that r 3 is the third compo- 
nent of total isospin and S is the (total) nuclear-angular-momentum operator. The 
first and second terms in Eqn. (20b) determine 3 He and 3 H, respectively, while upper 
and lower components refer to ± in Eqn. (20a). Finally we obtain 

4l = ^<^) + — W(2)-0, (22a) 

and 

*S = — (OS + - (r>S , (226) 

where the point-neutron, SU(4) symmetry limit (i.e., the S-state only, which is indi- 
cated by the arrow) is discussed below. 

This result is very easy to interpret. In the dominant S-state (corresponding to 
7 = and £ = 1) the two "like" nucleons (e.g., the protons in 3 He) have opposite 
and cancelling spins, while the "unlike" nucleon (e.g., the neutron in 3 He) carries 
all of the spin and determines both the magnetic moment and the single-nucleon 
contribution to the hyperfine structure (if we ignore meson-exchange currents and 
the convection current). Ignoring those currents leads to = r )ii p + £/i n — > ji n 
and yt/ajj = £/i p + 7/i n — > fi p , where the arrow indicates the SU(4) limit. The 3 He 
single-nucleon contribution to hfs in this limit comes solely from the tiny neutron 
contribution, while the corresponding 3 H contribution is simply (r)^ and becomes 
identical to the free proton Zemach moment. 

The two-nucleon contributions are more complicated and are determined by cor- 
relation functions. In the deuteron these correlations must be between a neutron and 
a proton and are of the types: e p — fi n , e n — fi p , and e n — e p (convection current only). 
In the trinucleons there are additional types: e p — fi p , e n — /i n , and e p — e p (convection 
current only). As we noted earlier the correlations involving e n will be very small, 
and we have ignored the tiny e n — e n convection- current correlation in 3 H. 

The two-nucleon contributions contained in Eqns. (15) and (16) can be manipu- 
lated into simpler forms by shifting x by x\ and y by x'j, leading to 

^Low — X! — 7pij{ x ij){3&j ■ XijXij — (Jj)\ , (23a) 
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Cow = — — - E CiiMLT , (23b) 
16M/i N ^ 

where = x (7^ — ffj) and the three correlation functions are defined by 

djir) =Jd 3 xJ d 3 ye t (x)fi :i (y)\x-y + r\ , (24a) 

= Pii^pPj + fi n hj)C DD (r) + Uiifippj + hJi^CdgIt) , (246) 

Cij(r) = J d 3 x J d 3 yei(x)ej(y)\x- y+ r| 3 , (24c) 

= PiPjC DD (r) + (ftihj + hiPj)C DG (r) , (24d) 
1 / d 2 1 d \ r r 
Ci3 ^ = 3 (rf^ ~~ rdr) J ^ X I ^V^^v)^ ~V+ ? l 3 ' ( 24e ) 

= Pi(/^ p Pj + n n hj)C DD {r) + hiiVpPj + /VbO^WO , (24/) 

where in Eqns. (24b, d,f) we have decomposed the charge and magnetic distributions 
in terms of isospin projectors and radial functions CddM and Cdg^)- Explicit forms 
for these functions are given in Appendix A. Note that the quantity in parenthesis in 
Eqns. (24b, d) is the (dominant) spin part of the magnetic moment operator, and that 
we ignore the contribution of two neutron charge distributions in Eqn. (24d). In the 
limit of no finite size each of the three non- vanishing radial functions (Cdd,Cdd,Cdd) 
simply equals r. 

The special case of the deuteron is easily dealt with. With our assumptions about 
the nucleon form factors there are only two distinct types of products contained in 
iiij and eifij] these are the dipole-dipole form of e p — ji n correlation (Cdd) and the 
dipole-Galster form of e n — e p and e n — correlation (Cdg)- The conventional form 
of the deuteron wave function (suppressing the spin and isospin wave functions) is 

^.(.W + I^W)^), (25a) 

which leads to the useful relations 

(4vrr 2 )^(r)i(o ; 1 + a 2 )^ d (r) = SA(r) + B(r) (ZrS ■ r - S) , (256) 

-SD(r) - 
2 v ; 4 



(4vrr 2 )^(r)L^(r) = -SD(r) - -D(r) (ZrS • t - S) , (25c) 



using L12 = 2L for the deuteron, while 

A(r) = u 2 (r) - \w 2 {r) , (26a) 

B(r) = -=u(r)w(r) + \w 2 (r) , (266) 
v2 

D(r) = w 2 (r) . (26c) 
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This leads immediately to 



4 2) = rdrA(r) (^C DD (r) + ^C DG (r)) 




^ ( ^C DD (r) + ^^(r)) + —D(r)C DG (r) , (27) 

where we have removed a factor of 1/2M from the magnetic moments (i.e., they are 
now given in units of nuclear magnetons). In the limit of vanishing neutron charge 
distribution and point protons, this expression becomes 

«F-y irr(Mr)-^f), (28) 

which is Low's expression [13] for the complete deuteron finite-size effect (in leading 
order). In the next Section we will refer to the integrals of radial functions such as 
A{t)Cdd(j'), (i-e., including the numerical factors, but without the magnetic moments 
in Eqn. (27)) as "Low moments." The two Low moments for point-like nucleons are, 
therefore, / rA(r) and — | / rB(r), as given in Eqn. (28). 



4 Numerical Evaluation 



The proton hfs has been discussed in detail recently[4, 15] and we have nothing more 
to add. The recently evaluated proton Zemach moment is listed in Table II, and 
it leads to a —58.2(6) kHz contribution to the hydrogen hfs, which equals —41.0(5) 
ppm. When added to the usual QED and recoil corrections [3, 4, 15] there is a 3.2(5) 
ppm discrepancy with experiment, which can be attributed to hadronic polarization 
and (possible) additional recoil corrections. 

We begin our discussion of the deuteron with the single-nucleon contribution given 
by Eqn. (19). Table II lists the proton Zemach moment and the neutron moment 
determined by our choice of form factors. The neutron result is only 4% of the proton 
value in magnitude, and the opposite sign reflects the fact that the (overall neutral) 
neutron has negative charge at large distances that balances positive charge at short 
distances. Using the value of Pd — 5.67% (corresponding to the AV18 potential 
model[23]) in Eqn. (19) leads to the nucleon-moment deuterium hfs contributions 
listed on the right side of Table II. The proton result differs from that in hydrogen 
by the factor of (1 — \Pd) and the statistical factors for the deuterium and hydrogen 
hfs. 
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Table II. 



Table 2: Neutron and proton Zemach moments and their single-nucleon contribution 
to the deuterium Is hyperfine structure. 



Nucleon Zemach Moments 


Zemach moments 


Deuteron Nucleon-Moment hfs 


proton neutron 


proton neutron total 


1.086(12) -0.042 fm 


-40.0 -1.1 -41.1 kHz 



Table III. 

Table 3: Deuterium Low moments from various parts of the nuclear current. The 
A-terms are the space-scalar contribution (first term in Eqn. (23a)) from the spin- 
magnetization current, the B-terms are the corresponding space-tensor terms (second 
term in Eqn. (23a)), while the D-term arises from the convection current (Eqn. (23b)). 



Deuteron Low Moments 


A p t B pt 


Add A D g B D d B dg D D g 


3.081 -0.115 


3.271 -0.015 -0.126 0.001 -0.003 fm 



The deuterium Low moments are listed in Table III, and the resulting hfs is listed 
in Table IV, both for point-like nucleons (only the proton charge contributes) and for 
nucleons with finite size. The moments themselves and the resulting hfs are defined 
in Eqn. (28) for point-like nucleons and in Eqn. (27) for finite nucleons. It is obvious 
that d-waves and the neutron's charge distribution play a minor role. The proton 
charge distribution and the neutron magnetic distribution have a somewhat larger 
effect; Add is larger than A pt by about 6%. 

One can also compute the Zemach moment of the entire deuteron by construct- 
ing the charge (F c h) and magnetic (-F m ag) form factors and using the equivalent 
momentum-space version of the Zemach moment formula: 

(r> (2 ) = ~ /f ^(i^g 2 )iWg 2 ) - l) • (29) 

Various contributions and limits are listed in Table V. Results for point-like nucleons 
are listed to the left and include the contributions from the s-wave spin-magnetization 
current, the d-wave spin-magnetization current, and the orbital (convection) current, 
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Table IV. 



Table 4: Contributions to the deuterium Is hyperfine structure from the Low moments 
compiled in Table III. 



Deuteron Low-Moment hfs 




Add A D g B D d B D g D D g 


84.9 -3.2 


90.2 0.6 -3.5 0.0 0.0 kHz 



followed by the total contribution. Including identical dipole nucleon form factors for 
the proton's charge and the neutron's magnetization densities (which multiplies both 
F c h{<l 2 ) and F mag (q 2 ) by F D (q 2 ) - see Eqn. (A3)) leads to the rightmost result and 
corresponds to an increase of about 10%. The experimental result of 2.593(16) fm 
was obtained directly from the electron-scattering data[15], and is approximately 2% 
smaller (4 standard deviations) than our non-relativistic calculation. This difference 
is the expected size of relativistic corrections from MEC. 

Table V. 

Table 5: Deuterium Zemach moments from various parts of the nuclear current. 



Deuteron Zemach Moments 


Point N 


Finite N 


Experiment 


L = L = 2 Orb 


Zemach 


Zemach 


Zemach 


2.324 -0.035 0.094 


2.383 


2.656 


2.593(16) fm 



The one-body (nucleon Zemach) and two-body (Low) contributions to the total 
deuterium hfs are listed in Table VI. Because there is then no point-nucleon contribu- 
tion to the one-body part, the Low term is the sole contribution and leads to a very 
large result. The finite-nucleon case has considerable cancellation between the two, 
and totals only about half the size of the point-nucleon limit. One can also break the 
total result down into deuteron Zemach (elastic) terms plus polarization (inelastic) 
terms. This is indicated in Table VII. The polarizability term is more than twice the 
elastic (Zemach) term, and reflects how easily a weakly bound system can be excited 
compared to a system like the nucleon, which is difficult to excite. In the latter case 
the polarization term is only about 10% of the Zemach term. 

The physics of the nuclear correction to the deuterium hfs is straightforward and 
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Table VI. 



Table 6: Contributions to the deuterium hyperfine structure from one-nucleon and 
two-nucleon operators and their total. 



Deuteron hfs - Nucleon + Low Moments 


Point N 


Finite N 


Nucleon Low Total 


Nucleon Low Total 


0.0 81.8 81.8 


-41.1 87.3 46.2 kHz 



Table VII. 

Table 7: Contributions to the deuterium hyperfine structure from elastic (Zemach) 
and inelastic (polarization) intermediate states. 



Deuteron hfs - Zemach + Polarization 


Point N 


Finite N 


Zemach Polar Total 


Zemach Polar Total 


-29.5 111.2 81.8 


-32.8 79.1 46.2 kHz 



completely dominated by the proton Zemach moment and the e p — fi n Low contri- 
bution. Nuclear structure plays little role except to fix the size of the (radial) Low 
moment. The signs were fixed by the sign of the proton magnetic moment (for the 
one-body term) and the neutron magnetic moment (for the two-body term), and are 
opposite. If we incorporate the additional minus sign in Eqn. (lib), the naively ex- 
pected sign of the one-body terms should be — , while that of the Low contribution 
should be +, as we found in the deuterium case. As we will see, however, nuclear 
structure can play an exceptional role in the trinucleon, and these expectations are 
not fulfilled in two cases. 

The required 3 H and 3 He matrix elements were calculated using wave functions 
obtained from a Faddeev calculation [22]. The (second-generation) AV18[23] potential 
was used, together with an additional TM' three-nucleon force [24] whose short-range 
cutoff parameter had been adjusted for each case to provide the correct binding 
energies. Individual one-body (labelled "nucleon") and two-body (labelled "Low") 
terms are tabulated together with their total in Table VIII. Note that the 3 He case 
(which has proton number Z = 2) is uniformly enhanced (compared to the H cases) 
by a factor of Z 3 = 8 contained in |0 n (O)| 2 in Eqn. (11a). For the same reason the 
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two protons in 3 He effectively double that Low moment. Taking those factors into 
account the 3 He Low term becomes comparable in size to that of the deuteron. 

The (approximate) SU(4) symmetry that dominates light nuclei[25] provides an 
explanation for the relative sizes of the entries in this table, as well as the unexpected 
signs (see above) of the 3 H one-body term and the 3 He two-body term. The two 
protons in 3 He have their spins anti-aligned in the SU(4) limit, and this cancellation 
leads to the small net result and unexpected sign for the one-body part, which is 
determined by small components of the wave function. The protons in X H and 3 H 
make comparable one-body contributions, since the proton in 3 H carries the entire 
spin in the SU(4) limit. For the same reason the Low term in 3 H is very small because 
the two neutron e p — ji n terms largely cancel, since the proton carries all of the spin 
in the SU(4) limit. 

The neutron Zemach moment plays only a very small role in the one-body terms 
(as we found for the deuteron) except for 3 He, which has a greatly suppressed proton 
contribution. The tensor term also is quite small. The convection current terms are 
negligible for 3 H, but the e p — e p contribution in 3 He is approximately 5% of the total. 

Table VIII. 

Table 8: Contributions to the trinucleon hyperfine structure from one-nucleon and 
two-nucleon operators. 



Trinucleon hfs 


3 H 


3 He 


Nucleon Low Total 


Nucleon Low Total 


-50.6 -9.6 -60.1 


14 1428 1442 kHz 



Our final results are listed in Table IX. The first line of the table is the same as 
that in Table I, showing the fractional difference of experiment and QED theory in 
ppm. That fractional difference is recomputed in the second line when the nuclear 
corrections are added to the theoretical result. In the proton case the Zemach and 
recoil corrections slightly over-correct, but the overall result is consistent with the 
expectation that the polarization corrections are positive and must be less than 4 
ppm[26, 27]. For the nuclear cases the quality of our results must be considered quite 
good, given the size of our hadronic expansion parameter. The deuterium case is 
particularly close to experiment, and this is likely due to the small binding energy, 
which tends to minimize relativistic corrections [6]. The quality of our trinucleon 
results range from very good in the 3 H case (~ 3% residue) to adequate in the 3 He 
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case (~ 20% residue). The large disparity in the two cases is undoubtedly due to 
missing MEC, particularly the isovector ones. Even this amount of missing strength 
is only slightly larger than our expansion parameter. 

Table IX. 

Table 9: Difference between hyperfine experiments and hyperfine calculations for the 
ls-state of light hydrogenic atoms, expressed as parts per million of the Fermi energy. 
The first line is the difference with respect to the QED calculations only, while the 
second line incorporates the hadronic corrections (Zemach moment for the proton 
and nuclear corrections for the nuclei) calculated above. 



(£ms P - E^)/E F (ppm) 



Theory 


H 


2 H 


3 H 


3 He+ 


QED only 


-33 


138 


-38 


-212 


QED + hadronic 


3.2(5) 


-3.1 


1.2 


-46 



Previous work on this topic is quite old[13, 14, 28, 29, 30, 31], except for the 
deuterium [32] case. The older work relied on the Breit approximation for the electron 
physics, which is sufficient only for the leading-order corrections. It used an adiabatic 
treatment of the nuclear physics based on the Bohr picture of the nuclear hyperfine 
anomaly, which is far more complex than the treatment that we have presented. 
Uncalculated QED corrections and poorly known fundamental constants (such as a) 
led to estimates of nuclear effects that were many tens of ppm in error. Although 
the nuclear physics at that time was not adequate to perform more than qualitative 
treatments of the trinucleons, the SU(4) mechanism was known and this allowed 
a qualitative understanding. The only previous attempt to treat the three nuclei 
simultaneously was in Ref. [31]. They found nuclear corrections of about 200 ppm 
for deuterium, 20 ppm for 3 H, and —175 ppm for 3 He + . Except for the deuterium 
case (which involves significant cancellations) this has to regarded as quite successful, 
given the knowledge available at that time. 



5 Conclusions 



We have performed a calculation of the nuclear part of the hfs for 2 H, 3 H, and 3 He + , 
based on an expansion parameter adopted from x?T, a unified nuclear model, and 
modern second-generation nuclear forces. This is the first such calculation, and the 
results are quite good. Details of the results can be understood in terms of the 
approximate SU(4) symmetry that dominates the structure of light nuclei. 
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6 Appendix A 



The correlation functions that we require are built from the charge and magnetic form 
factors for protons and neutrons. The generic correlation function has the form 

P(2)(z) = J d 3 x p ch (x) p mag (\x + z\), (Al) 

where the charge density p c h( x ) and magnetization density p mag (y) are normalized to 
one. Although a variety of functional forms have been proposed for these densities, 
few published forms have high accuracy over the low-momentum-transfer region that 
is important for Zemach moments. Fortunately the proton's Zemach moment was 
recently determined to high accuracy, and the neutron's is sufficiently small that any 
credible model should suffice. For the neutron charge form factor we will assume a 
modified Galster form [16] 

W) = ^, (A2) 

with (3 = 4.0285 fm _1 (determined below) and A = 0.0190 fm 2 . This is accurate 
enough[19] for our purposes, both at low values of q 2 and at moderate values of q 2 . 
For the proton charge form factor and the proton and neutron magnetic form factors 
we choose the tractable and venerable dipole form 



which is a reasonable (but only moderately accurate) approximation. 

These two forms inserted into Eqn. (Al) generate the two correlation functions 
that we require: Pdd( z ) and Pdg( z ) : 

Anp DD (z) = ^ exp (-/3z)(3 + 3((3z) + (f3z) 2 ) , (A4) 



and 

Air r>n^( z) = 

384 



^p DG (z) = ^ exp (-p z ) (9 + 9(0z) + 2{Pz) 2 - (f3zf) . (A5) 
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The first moment of these functions is the linear Zemach moment 

( r )?2) = J d 3 rrp DD (r) = j d 3 x J d 3 y p D {x) p D {y)\x - y\ = ^, (A6) 

which we identify with the recently determined proton moment: 1.086(12) fm. This 
restricts (3 to be 4.029(45) far 1 , which we also use for the neutron. The rms radius for 
a dipole with this value of f3 is 0.86 fm, slightly smaller than the proton's charge radius, 
but slightly larger than the magnetic radius by a few percent, and this represents our 
level of accuracy (except for the measured proton Zemach moment). 

We can use these functions to determine the appropriate correlation functions 

Cdd{t) = J d 3 x Pdd(x)\x+ r\ 

8 exp(-/3r) / 8 29 5(/3r) (/3r) 2 \ 
~ T + JjPr) \Jj3rj + ~8 + 8 + 24 J 

-| + "'- < A7 > 

where the limiting form holds only for small (f3r), and similarly 

C D G(r) = J d 3 xp DG {x)\x + r| 

= x > 3 (" W) + (m + 279 + 87(W + u(,3rr + ij3r)3 )) 

35X13 +-... (A8) 



64 

The remaining functions that we require are determined by 

1 3 



C DD( r ) = J d 3 xp DD (x)\x + r| 



a 1 / 240 in/n . , „ N /240 165 21(/3r) (pr) 2 \\ 

= r+ p{W) + 48(5r) " exp ( -' 3r) ( (M + - + ^ + 2 J J 



2/5 3 4/5 

and 



CdgM = y d 3 xp DG (x)\x + r\ 



I w - M + t K + 975 + 207(w + 22W2 + w3 ) ) 

945A 35A/5r 2 , A _ N 

1 . (A10) 

16/3 32 v ; 
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From the former we obtain 



r 



16 80 exp(-/3r) / ,„ . lln n2 19(/3r) 3 (/3r) , 

= r + Wr) ~ W + r + 8 ° ( " r) + 24( " r) + 6 + 6 j 

+ (AH) 

and 

/3r 2 

(A12) 

and from the latter 



= V (" W) + WY> ~ "W 2 ( 40 + 4 ° (,3r) + WW + « + 48 + 48 
-A/j(-| + ...), (A13) 

and 



= A/? 



3 \ c?r 2 r dr / 
4 120 



(/3r) {f3rf 



+ ( 120 + 12 °(^) + 56 (^) 2 + 16 (^) 3 + 3 (^) 4 + lz |r ! + ^r) 

^,(^! + -). CAU) 

Note that these functions have been normalized so that Cdd, Cdd, andCoD become 
r in the limit of large f3. 

We resort to the simple zero-range approximation in order to make a rough esti- 
mate of the effect of nucleon finite size on the dominant Low moment of the deuteron. 
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This approximation is most accurate for asymptotic (long-range) quantities, but will 
substantially overestimate short-range effects. We ignore the d-state and assume 
everywhere the asymptotic s-wave function (iV r e~ Kr / ' \f^nr where k ~ 0.235 fm _1 ). 
This leads to an expansion parameter x = 2k/ (3 ~ 0.115. The matrix element of 
Cdd (relative to the point-nucleon value) is approximately (1 — x 2 (y + 8 In (x) + . . .), 
which produces an increase ~ 18% from nucleon finite size. This is too large by a 
factor of two compared to detailed calculations, but shows that the finite-size effect 
is enhanced by the large numerical coefficient of the logarithmic term beyond what 
is expected from a 0(x 2 ) ~ l/(3 2 correction term. 



7 Appendix B 

Many features of the trinucleon systems, 3 He and 3 H, can be determined in a semi- 
quantitative fashion (accurate at the ~10% level) by simplifying the wave functions 
to the dominant component alone. Wave function components have traditionally 
been classified according to their combined spin-isospin symmetry, determined by the 
generators of SU(4). In this scheme the SU(4) generators for a system of A nucleons 
are determined by the intrinsic spins and isospins of the individual nucleons: 

A 



and 

yrka 



E°?> (Bl) 
i=\ 

X>°, (B2) 

X>>*. (B3) 



i=i 

All wave function components are labelled by the (combined) intrinsic spin of the three 
nucleons (S = 1/2 or 3/2) and the total isospin (T = 1/2 or 3/2). Wave function spin 
and isospin components are then determined by: (1) the [4] or antisymmetric state 
(S =1/2, T—l/2), which combines with a completely symmetric space wave function 
to form the dominant S-state; (2) the mixed-symmetry state, which can be separated 
into (S = 1/2, T = 1/2), (S = 3/2, T=l/2), and (S = 1/2, T = 3/2) components. The 
first term contributes to the small S'-state, while the second generates the D-state(s), 
and the third contributes only to tiny isospin impurities. The remaining spin-isospin 
representation has a tiny (S = 1/2, T = 1/2) symmetric component (called the S"- 
state, with a totally antisymmetric space wave function), and a (S = 3/2, T = 3/2) 
D-wave isospin impurity. The number of components in the order discussed is (4) + 
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(4 + 8 + 8) + (4 + 16) = 64, as expected. Ignoring the tiny S"-state, very small 
isospin impurities, and the negligible P-states, we can therefore write in an obvious 
but schematic notation for the trinucleon wave functions 

tttri = S © S' © D . (B4) 

It was shown many years ago [20] that expectation values of the SU(4) generators for 
the trinucleons have very simple forms. 

(J2a l )=2S(l-2P D ), (B5) 
i=i 

and 

(t *0 = "2^3 (l - - \P D ) , (B6) 

where a spin-isospin expectation value of the nuclear spin and isospin operators is 
still required. 

As expected, the spatially symmetric S-state dominates the trinucleon ground 
states (Ps ~ 90%) because it minimizes the kinetic energy. The mixed- symmetry 
S'-state is much smaller (Pg/ ~ 1%), while the very strong nuclear tensor force gener- 
ates a relatively large D-state component (Pd ~ 9%). If one ignores the S'-, S"-, P-, 
and D-state components, the remaining S-state wave function factorizes into a com- 
pletely symmetric space wave function and a completely antisymmetric spin-isospin 
wave function, which greatly facilitates calculating matrix elements. The mixed spin- 
isospin generator Y ka has the very useful and simple property for the S-state 

Y ka \S) = -2S k r a \S), (B7) 

which follows (except for the factor of (—1)) from the Wigner-Eckart theorem and 
the properties of the [4] state. 
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